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Abstract 

The system of hydrodynamic-type equations, derived by two-side distribution 
function for a stratified gas in gravity field is applied to a problem of ultrasound 
propagation and attenuation. The background state and linearized version of the 
obtained system is studied and compared with the Navier-Stokes one at arbitrary 
Knudsen numbers. The WKB solutions for ultrasound in a stratified medium are 
constructed in explicit form. The problem of a generation by a moving plane in a 
rarefied gas is explored and used as a test while compared with experiment. 

1 Introduction 

Recently the problems of Kn regime wave propagation was revisited in connection with 
general fluid mechanics and nonsingular perturbation method development PjJ 12 El H] • A 
generalized Boltzman theories Ej also contributed in a progress with respect to this 
important problem. 

In |Zj the propagation of one-dimension disturbance was studied on the base of the 
method of a piecewise continuous distribution function launched in a pioneering paper of 
Lees jH] and applied for a gas in gravity field in [TJ We derived hydrodynamic-type 
equations for a gas perturbations in gravity field so that the Knudsen number depends 
on the (vertical) coordinate. The generalization to three dimensions is given at [12J 

The derivation of the hydrodynamic-type equations is based on kinetic equation with 
the model integral of collisions in BGK (Bhatnagar - Gross - Krook ) form which collision 
term is modelled as v (fi — /) , via local-equilibrium distribution function /; and the 
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non-equilibrium one is expressed as / + at v z > 0, and as / at v z < 

, (V-U*) 2 . 



f ~ 3/9 ±3 eX P(" -!-2 



the v T = \j2kT/m denotes the average thermal velocity of particles of gas, v = v(z) is 
the effective frequency of collisions between particles of gas at height z. It is supposed, 
that density of gas n, its average speed U = (u x , u y , u z ) and temperature T are functions 
of time and coordinates enter the local-equilibrium fa. The resulting system is 
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where 

777 — * — * 777 

.h = -<(y z -uy{v-uy>, j 2 = -<(v z -uy> . ( 2 ) 

The increase of the number of parameters of distribution function results in that the 
distribution function differs from a local-equilibrium one and describes deviations from 
hydrodynamical regime. In the range of small Knudsen numbers I << L we automatically 
have n + — n~ , U + = U~ , T + = T~ and distribution function reproduces the hydrodynam- 
ics of Euler and at the small difference of the functional "up" and "down" parameters 
- the Navier-Stokes equations. In the range of big Knudsen numbers the theory gives 
solutions of collisionless problems [Hj. 

We used a set of linearly independent eigen functions of the linearized Boltzmann 
operator, that in the case of the BGK equation is: 

<fi = m , (f4 = m(V z - U z ) 2 , 

tp 2 = mV z , (f 5 = -m(V z - U Z )\V - U\ 2 , ( 3 ) 

<^3 = ^m\V-U\ 2 , = ^m(V z -U z ) 3 . 
Let's define a scalar product in velocity space: 

< tp n , f >=< ip n >= [ dv <p n (t, z, V)f{t, z, V) . (4) 
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< ip x >= p(t, z) , < V2 >= pU , < y? 3 >= , 

< V?4 >= P zz , < </? 5 >= ?z, < V?6 >= Qz- 

Here p is mass density, P 2Z is the diagonal component of the pressure tensor, q z is a 
vertical component of a heat flow, q z is a parameter having dimension of the heat flow. 

The system Q of the equations according to the derivation scheme is valid at all 
frequencies of collisions and within the limits of the high frequencies should transform to 
the hydrodynamic equations. 

If we estimate the functions j-^ as small, that corresponds to small Mach numbers 

U 

M = max\ — I . We shall base here on an expansion in M, up to the first order. In this 

Vt 

approach the functional parameters of the two-fold distribution function 

n + = n (l + anf) , n~ = n Q (l + an±) , p = nm 
V+ = V (l + aV 1 + ) , V- = V (l + aVf) , 
U+ = aVvUt , U~ = aV U{ 

Let's evaluate the integrals (j2j) and © directly, plugging the two-side distribution 
function. In the first order by Mach number a 



n = n + \-j=~ + 2 n o n i + 2 n o n i ) a 

U = (-|^ - ^ + \VJJt + ^ + lV Ur + \ V ^f) a 

l^pT = \n V * + (fno^o 2 ^ + fn K, 2 < + §n V 2 nr + n -^f- + f^o^f - n -^B-) c 
iP zz = |n y 2 + (jnoVo V - + §fio*W + JnoVoV + faV *VT + =2*J£) a 

1„ _/ / 5„T / 3rr+ 5 „ y 2tj , 5„ T/ 3rr- 3 noVb 3 V~ 1 noV 3 n+ 1 «o%V , 3 ^0^+ \ 
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(6) 



If — 5 r,T/ 4 4-/' 5 « 1/4^+ 1 5„ t/4t/- 3 "oVb 4 ^ , 5^ T/4 T/ + , 5 T/ 4 - , 3 "oK> 4 ' 

-Ji - ^n V + I Y6 n i + l n o y o v x - 5 — ^ h 4^0 ^0 V^i + 16^0 ^0 % + 2 — 7^ 

i h = + l±n V ±ni + |n ^Vf + > ^oV - + f + 

(7) 

Solving the system (JHJ), we obtain for the parameters of the two- fold distribution 
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function: 

+ _ 3 Uy/n n P zz kpT q zy /ii q zy /ii | 3 nU^/n 

2 Vo n mnoVo n m 2 Vo mnoVo mnoVo 2 n Vo 

= 3 Uy/7r + 3 ^ + 3 fc P T + 7 JzVn_ _ 3 jW^L _ 3 nU V^ 
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^ + _ 1 raE/y? _ 1 3 fcpT _ l_n q z ^n lUyft _ g z ^n 
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£/ 3 \fnP zz 3 \fnkpT q z q z 

Vo 2 mnoVo 2 riom 2 Vo mra Vo mnoVo 

The values of integrals (J2J) as functions of thermodynamic parameters of the system 
(0) are linked to the thermodynamic variables as: 

5 ,kT . 9 llkT P zz 9 . k . 9 mm 

Jl = ~2 P{ } + + 4 ( m } PT ° T ' 

3 # 2 9^ 3 f 2 ( 8 ) 

h = -~P{ ) + 7 + t(— ) P?oT. 

2 m 4 m 4 m 

So we have closed the system hence a modification of the procedure for deriving 
fluid mechanics (hydrodynamic-type) equations from the kinetic theory is proposed, it 
generalizes the Navier-Stokes at arbitrary density (Knudsen numbers). 

Our method gives a reasonable agreement with the experimental data in the case of 
homogeneous gas [7j. In the paper [7j the expressions for (Jip) are obtained with account 
some nonlinear terms, that finally lead to more exact results. 

2 Stationary case (undisturbed atmosphere). 

Let's linearize the system (JTJ) this way: 

p = poO)(l + epi(t, z)) , U z = ell z i(t, z) , 

P(t,z) = P (z)(l + eP 1 (t,z)) , T = T (z)(l + sT 1 (t,z)) , 

q z = q z0 (z)(l + eq zl (t,z)) , q z = q z0 (z)(l + eq zl (t, z)) , e«l. 
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We obtain in the zero order: 
PoU) 

Tz qo{z) = ' 

a£ m 

5 fc m , . d i , -Pnl- 2 ) d ^ . . 11 A; „ . . c? _ . . 

+~_T * _P z - ^H^PoW + T -Po * h-To * = 0, 
4 m 02; Po(-2J u- 2 4 m dz 

9 k d 3 k 2 d 

u(z)q z0 (z) + -—P (z)—T (z) - -—p (z)T (z)—T (z)- 

3k2 T 2 (z) dm o(z) + ^ k T(z) d 7(z) 3 P °(*>V°(*> -Q 
l/» ( ^ P ° W + 4m To{z) Tz P ° {z) ~ 2—JJz) " °- 

Some version of such system that leads to a non-exponential density dependence on height 
was studied in [TUJ , the paradox was discussed at jTT] . 
Let's solve the zero order system. 

q zQ = C X . If P = ^p T , then q z0 = C 2 = \C X . If C x = 0, then T = C 3 = const and 
we'll have exponential density dependence on height. 
We obtain in the first order: 
d d V T 

m» +VT a; Ui -i Ui = > 

dU x V T dP x V T( p , n 

dpx dTx 5 dU x 2 dq x V T 2 
^ + ^ + 3^ + 3^-^ (f/l + 3 gi) = °' 

-P X + 3V T ^ + 2V T -^ - -^(2q x + Ux) + u(Px - px - Tx) = , 

dq x 1 dP x 9 dT x 1 dp x 3V T 

5F + 8^ + 8^-8^-8¥ (Pl + Tl - Pl) + Z/gi = ' 
dq x 3 dP x dT x dp x 3V T _ 
^ + 8 M ^z- + -dz~- a7 ) + 8F (Pl - pi - Tl) + Z/9l = °- 

3 Construction of solutions of the fluid dynamics sys- 
tem by WKB method. 

In this section we apply the method WKB to the system (|l()jl . We shall assume, that 

on the bottom boundary at z = a wave with characteristic frequency uj is generated. 

Next we choose the frequency ujq to be large enough, to put characteristic parameter 
3uj H 

£ = ^> 1. We shall search for the solution in the form: 

v T 

M n = ^ n exp(ia;ot) + c.c. , (11) 



(10) 
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where, for example, tpi, corresponding to the moment Ml, is given by the expansion: 

6 oo i 

^i = EE 7^ A t ] , (12) 

fc=im=i y^) 

here (fk(z) - the phase functions corresponding to different roots of dispersion relation. For 
other moments M n , n — 2, ... ,6 corresponding functions ip n are given by similar to (|T2"j) 
expansion. The appropriate coefficients of the series we shall designate by corresponding 
C$ D$ E$ F$ . Substituting the series (O at the system (JTUJ) one arrives at 
algebraic equations for the coefficients of (J12)) in each order. The condition of solutions 
existence results in the mentioned dispersion relation: 



(13) 



&rf> + (-fin - | + §m 2 ) r? + (-in 3 + fin - f m 2 + f ) rj- 

-1 - 3iu + 3u 2 + iu 3 = 
Here for convenience the following designations are introduced: 

{-dF) =15^' M = ^ exp( "^' 

where z — — . For the coefficients . . . the algebraic relations are obtained: 

H 

{k) VSO A[ k) {k) _ 1 4 fc) (-25 + 20m + 37/) _ 5 4 fc) 

1 5 ^ 1 ) 1 5 

6 3 —10 + lOzii + 97/ 3 

(fc) 5 AjVg(l + 3^) 

1 12 ^(-10 + 10m + 9t/) ' 

(fc) , 1 /30 (50 - 100m-135r]-50M 2 + 190m?7 + 81r7 2 + 50M 2 r/-30^ 2 u)^ {fc) 

1 36VV ^(-10 + 10m + 9t/) 1 

The dispersion relation ()13j) represents the cubic equation with variable coefficients, 
therefore the exact analytical solution by formula Cardano looks very bulky and inconve- 
nient for analysis. We study the behavior of solutions at v — > (free molecular regime) 
and v —>■ oo (a hydrodynamical regime). 

At the limit of collisionless gas v — the dispersion relation becomes: 

54 „ 63 2 . 1* 



-77 V H H -1 = 0. 

125 ' 25 ' 5 ' 
The roots are: 

7/! « 3.80 , 7/2 « 0.37 , 7/3 « 1.67 . 

In a limit z/ — ► oo (a hydrodynamical limit) for specifying roots (fT3*|) by the theory of 
perturbations up to u 3 for the three solutions branches it is obtained: 

m = 1.00 - 2.32u- 2 + ^l^Oir 1 - 4.88u- 3 ) , 

t/2 = 71.67m + 2.33 + 0.64m" 2 , 

r/ 3 = -1.39m 2 + 2.50 + i(3.89u - 1.20m- 1 ) . 
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The first root relates to the acoustic branch. Accordingly, for the ki t ± = ±yjr]i we have: 



k 2 ,+ ~ 



1.00 - 0.98m~ 2 + ^O.eOn- 1 - 1.85u- 3 ) , 

v/w(l - ^(O.G^- 1 + 0.019m- 3 ) + v/w(l + z)(0.91 + 0.22«- 2 ) 

1.65 - 0.64m" 2 + z(1.18w + 0.094m- 1 ) . 



The solution of the equation (|T^|) at any u is evaluated numerically. As an illustra- 
tion let us consider a problem of generation and propagation of a gas disturbance, by a 
plane oscillating with a given frequency ujq. We restrict ourselves by the case of homoge- 
neous gas, because it is the only case of existing experimental realization. We evaluate 
numerically the propagation velocity and attenuation factor of a linear sound. 
1 
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Fig. 1. The inverse non-dimensional phase velocity as a function of the inverse Knudsen 
number. The results of this paper-1 are compared to Navier-Stokes, previous our work 
[7J-2 and the experimental data of Meyer- Sessler |T3j-circle. 
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Fig. 2. The attenuation factor of the linear disturbance as a function of the inverse 

Knudsen number. 



4 Conclusion 

In this paper we propose a one-dimensional theory of linear disturbances in a gas, stratified 
in gravity field, hence propagating through regions with crucially different Kn numbers. 
The regime of the propagation dramatically changes from a typically hydrodynamic to 
the free-molecular one. We also studied three-dimensional case [12] • Generally the theory 
is based on Gross- Jackson kinetic equation, which solution is built by means of locally 
equilibrium distribution function with different local parameters for molecules moving 
" up" and " down" . Equations for six moments yields in the closed fluid mechanics system. 
For the important generalizations of the foundation of such theory see the recent review 
of Alexeev [5]. 
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